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Abstract 

Hard-core final state interaction effects in the response function are investi- 
gated in the asymptotic limit of momentum transfer g — > oo. The relevant 
scattering contribution is displayed and a modification of the otherwise ex- 
pected free response is obtained for Bose and Fermi systems. A comparison 
with other treatments of final state interactions is made. 
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Important information on Bose and Fermi systems can be obtained from inclusive scat- 
tering of a weak probe in which energy uj and momentum q are transferred to these systems. 
Such measurements in the large q region offer the possibility to obtain the single-particle 
momentum distribution n{p). Indeed, if for large q the response function S{q, uj) approaches 
the impulse approximation (lA), then the momentum distribution can be readily obtained. 
In terms of the reduced response 0(5', y) = {q/m)S{q,uj) the lA is given as 



where y = muj/q — q/2 is a scaling variable, A the number of particles in the system, m 
the mass of the particle, and pi is the one-particle density matrix in coordinate space. In 
particular, an examination of Eq. (|ID implies that the lA gives a distinctive delta function 
peak for the case of a momentum distribution containing a condensate fraction no, namely 
n{p) = rincij)) +'^c(p), where ndj)) = no{2n)^6{p). However, final state interactions (FSI) 
can modify Eq. (|1|) and it is of special interest to study their effect on the condensate peak. 

The treatment of FSI is a very complicated problem in general and it is desirable to 
study the response in the asymptotic limit q —>■ oo, and y=const, where the problem is 
greatly simplified. Indeed, Eq. represents the lowest order term in a multiple scattering 
series (MSS), i. e., (p{q,j) — (p^^iy) = 0{v/q), where v is the potential. This means that 
the lA is exact in the asymptotic limit in the case of regular potentials. In addition, if the 
interactions contain a strong repulsive component (hard-core), only this part has to be taken 
into account in the asymptotic limit, while the regular part will not contribute. 

In the following we treat the asymptotic response in terms of a MSS describing the FSI. 
Being able to sum over the MSS into simple analytical expressions, we investigate for both 
Bose and Fermi systems how the lA is modified in this limit. In particular, we explicitly 
show how the 5-function condensate peak is smeared due to FSI. We also make a comparison 
with other approaches for the response function in the asymptotic limit, thus clarifying in 
this way the nature of the approximations underlying these. 

The exact asymptotic expression for the Green's function entering the response in the 



0^^(y) = / 
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case of hard-core FSI was obtained in Ref. [T[] by use of the well-known result that the 
Schroedinger equation coincides with the equations of geometrical optics when the energy 
of the struck particle goes to infinity. The geometrical optics solution is known for a wave 
function describing hard-core scattering of the knock-on particle '1' off spectator particles, 
'2',...,'y4'. It is a plane wave describing particle '1', modified by a subtraction of "shadow" 
regions behind the scatterers and adding the refiected waves |@]. The resulting form for the 
response 0(y) = limg^oo0(Q'7 y) is 

kj) = -^^e I d%,dz,dz[d\2...d'rAPAe{z[ - zi)e^y(^i " ^i) n(l - Su), (2) 

where 

Sii = e{a-\bi~ hi\)e{z[ - Zi + Wi)9{wi - zi + Zi) (3) 

with Wi = (a^ — \bi — bj^)^/^, Vi = (6i, zi), r[ = (bi, z[), and a is the radius of the core. pA 
is the y4-body density matrix defined as 

Pa = pA{ri, ra, va] r[, ra, va) = Al'^oiri, ra, rA)%{r[, ra, va). (4) 

where \E'o is the ground state wave function of the A-particle system, normalized in a box of 
volume V as (\E'o|^E'o) = 1- The corresponding n-body density matrices are defined by: 

Pn{ri, r„; r[, r„) = ^ J d^r^+i ■ ■ ■ d^rAPA^ri, r^; r/, r^). (5) 

Notice that the pa needed is off-diagonal in the coordinate of the struck particle Vi, and 
diagonal in the coordinates of the spectator particles. 

With its rather simple appearance, Eq. (H) gives the exact account for the FSI of the 
knock-on particle with spectators in the asymptotic limit. It is only the hard-core component 
of FSI which survives in this limit and modifies the free eikonal Green's function in Eq. (0) 
by the product of (1 — ^ij). This product forbids propagation along the hard-core shadow 
configurations (the hatched region in Fig. 1). By expanding the product in the integrand 

(I) 



Y\_{i — Sli) — 1 — ^^Sii + ^^SiiSij — ^ SiiSijSik + 

i i>j i>j>k 



(6) 
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one gets parallelly for the response function 



0(y) = E<^^(y) 



(7) 
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where the first term, 0i(y), corresponds to lA, Eq. (|T]) and the higher order terms represent 
a Glauber-type MSS, describing FSI of the knock-out particle with A — 1 spectator particles. 

As a first instructive example for the application of Eq. (|^) we consider the asymptotic 
response of a system of uqA particles occupying the zero momentum state, as can be obtained 
from the Bose condensation phenomenon in an ideal gas. We take the A-hodj density-matrix 
in the form = PaI = noAl/V"^, demanding pidri, r[) = n^p, where p = A/V . By substi- 
tuting p^] into Eq. (|^), we notice that its independence of the particles' coordinates enables 
Eq. (H) to be written as products of equivalent integrals over the spectator coordinates: 



^i^Hy) = / dzidz[(fbip9{z[ - zi) exp [iy{z[ - zi)] 



l--pE{z[-z,: 



A-l 
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with 



E{z[ - zi) = / d%idzi9{a - \bi - bi\)6{z[ - Zi + Wi)6{wi - zi + Zi) 



ra I 4a 

271 / dbiibii{z[ -zi + 2Ja^ - bj,) = na\z[ - zi + — ), 
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0(1) (y) = !^Re / dxe'y"^ 

TT JO 



1 2/ 4a, ^^-1 
1 - -pvra [x+Y^ 



(9) 



and bii = \bi — bi\. The function E depends only on x = z[ — zi and then Eq. (H) may be 
rewritten as 



(10) 



By expanding the integrand in Eq. (|10D using the binomial formula one explicitly reproduces 
the lA term and the MSS, Eq. (^. In the limit A ^ oo Eq. (p!OD goes over to 



0r^''(y) = ~6 3 Re / dxexplijx — pax] = —e ^ - — —j^ — -, 
71 Jo 7X {<^P) + y 



(y) = — ^ ' 

where a = vra^. 



(11) 



This result explicitly shows smearing of the (5-function peak due to FSI, first suggested 
by Hohenberg and Platzman 0. Notice that the 5-function peak survives in the lA for finite 
A, and disappears only if A ^ oo. This explains why smearing of the condensate did not 
occur in a truncated MSS |Q. Our treatment implies that it is not possible to determine a 
finite number of collisions leading to the extinction of the singularity as imprecisely stated 
in Ref. In addition, unlike Ref. P], in our treatment the role of the hard-core in the 
smearing is specific to the asymptotic limit. However, at finite q regular interactions can 
also produce the effect 



Eqs. (p!0|) , ([Tl|) give the asymptotic response for systems with a constant pA in space. 
Actually, such cannot be produced because of the presence of configurations with parti- 
cles' cores interceptions (square-hatched region in Fig. 1). In fact, the resulting Lorentzian 
response in Eq. (|l2) violates sum rules, and in particular, / 0^^^(y)c?y 7^ no, and in fact, the 
latter integral and / 0c^^(y)y^'iy diverge. These problems are overcome by modifying the 
density-matrix p^^, imposing the hard-core constraint. Specifically, we consider the A-body 
density-matrix 

P^Il=P^llvCir,,r'^), (12) 

where C(rj,r^) = Ylf=2^i''^ii ~ ~ ru = |rj — ri|, r'u = |r.j — r[\. Here the 

normalization factor rj is included, in order to to satisfy the condition for the one-body 
density pic = uqp, Eq. (||). We obtain 

v{x) = il-uix)/Vf-^, (13) 



^ Since the main contribution of the integral in Eq. (^) comes from any interval [0,x] with 
X >> l/{pa), the response could be approximated by contributions from a finite number of collision 



terms N ~ pax (replacing e^P^^ in the integrand in Eq. (11) by its Taylor expansion to A^th order). 
However, contributions from the interval [x, 00] remain singular unless one takes N ^ 00 which 
corresponds to considering an infinite number of collisions. 
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where 

8 14 
u{x) = -na^ + {na^x - ^^ttx^ - -TTa^)9{2a - x). (14) 

Configurations with particles' core interception for the hit and spectator particles, (Fig. 1), 
are actually excluded automatically by the (1 — Su) factors in Eq. (|^). Explicitly, when the 
latter Eq. is evaluated by substituting into it p^^l given in Eq. (0), the factor C(rj,r^) is 
redundant. For this density-matrix it immediately follows 



c 



(y) = — e 3 Re / dxexplijx — p{ax — u{x))] (15) 

TT JO 



The behavior of being other than Lorentzian at large |y|, allows for the (partly 

asymptotic) sum rules to be satisfied. Indeed one obtains from Eq. (P^) 



dj4>^'\y)=no (16) 

dY4>^'\y)Y = (17) 

Jdy4>^'\y)f=^-^<K> (18) 

where 



'2 



<K>=-^fj^\0,x) 
2m dx"^ 



(19) 



is the average kinetic energy. 

At low densities or at low momenta |ya| < 1 the second term in Eq. ([I^ ) can be neglected. 
In the limit A —>■ oo, rj ~ e 3 , and this Eq. gives 

7(2)/ \ ^0 Iz^Kif. ap 

0r(y) - — e 3 . 2 . (20 

TT (orp) +y 

which differs from the previous result in Eq. ( pA]) by the sign in the exponent. We notice also 
the higher peak modification as compared with the Lorentzian form that can be obtained 
from Ref. ||^ for a structureless fluid. 

The above simple model for the response reproduces asymptotic sum rules and embodies 
the smearing effect of the hard-core interaction by an infinite number of spectators. Other 



models with more realistic density matrices, are now shown to share yet the same basic 
underlying structure and behavior. 

The wave function in the product form \E'o = Ylij^j f{^ij) /'^^^'^ describes more suitably 
the correlations between the particles and is often used in calculations, (see for instance 0). 
The function f{rij) is directly related to the pair-correlation function and tends to 1 for 
larger than the correlation length. Correspondingly, the off-diagonal density matrix can be 
written as 

pt^ = ^X{f{r^^)f{r'u) J] f{r,,) (21) 

Now we substitute Eq. (pTj) into Eq. (|^) and use the expansion (^, (|^). The first term, 0i 
is the lA, Eq. (|l]). For the second term one finds 

0^*"^)(y) = — L-Re r dx [ (fri exp{ijx)E{ri, ri + xz) (22) 

T[A Jo J 

where z = q/q, x = z[ — Zi, and 
E(ri, r[) = J dh2pt\ri. r^, r[, r2)e{a - \b, - b2\)e{z[ - + W2)e{w2 - z, + z^). (23) 



Here p2 is two-body density matrix, Eq. (^. Notice that the 6'-functions in Eq. (|23|) reduce 
the r2-integration to the hatched region shown in Fig. 1. 

For the third order term, 03 also the integration over the spectator coordinates r2, Ts, is 
restricted by the same hatched region (Fig. 1). While the final expression now includes the 
three-body density matrix it can yet be approximately reduced to one involving two-body 
density matrices by noticing that the function /(r23), describing two-body correlations in 
Eq. (l2l|) , is close to one except for the region of small r23. For x large enough (larger than 
the correlation length), r23 will also be large for most values of the spectator coordinates, 
r2, r^, in their non- vanishing integration region (Fig. 1), and one can replace f{r2z) by 1 in 
the corresponding integral for 03. One obtains from Eqs. (^), (^T]) 



Using this result one gets for (/)3 

1 



(y) 



Re / dx d riexp(iyx) 



{A-2)E^{ri,ri + xz) 



2pf"'\ri, ri + xz) 



(25) 



vrA Jo 

Applying the same approximation for higher order terms we finally obtain for the asymp 
totic response 



0(^'^)(y) 



— -Re / dx I d^ripi^\ri;ri + xz) ex.'piiyx) 
TxA Jo 



1 - 



E{ri, ri + xz] 



-I A-l 



{A-l)p^f''\ri-n + xz)_ 
Here too, similarly to Eq. (|), the FSI modifies the lA result, Eq. (P, by an additional 
correction factor in the integrand. If pi{ri; ri +xz) const for x oo, the FSI correction 
factor would smear the 5-function peak in the response in the limit of A oo, just in the 



same way as in the previous example. Notice that the integrand (pq ) is exact for x — > cxo, 
and therefore the smearing would always occur in the asymptotic limit. 

Other treatments of the response can be compared with our asymptotic results. In 1973 
Gersch and Rodriguez (GR) obtained the following formula for the response, using the 
first cumulant term: 

P2{r,0;r + xz,0) 



1 /"OO P 

0^ (y) = — Re / dxpi{0;x) exp iyx + / d'^r- 
Tip Jo J 



pi{0;x) 



-ni<^^)(^r,r + xz) 



(27) 



where 



r, r + XZ] 



exp 



— / dx'[v{r + x'z) — v{r + xz)] 
q Jo 



(28) 



and V is the inter-particle potential. Here translational invariance is assumed, so that 
Pi{ri;r[) = pi(0;x), and P2{ri,r2;r[,r2) = p2{r,0;r + xz,0), where r = Ti - r2, and 
r' = r[ - r2. 

Although the original GR treatment was not intended for singular interactions |]T0| , Eq. 
(p7|) can be formally applied in the case of a hard-core interaction by taking w — > oo when 
the particles are within their hard-core range. We consider Eq. (^) in the limit q ^ oo 
and we take into account only the hard-core component of v in Eq. (^) (since the regular 
component does not contribute in this limit). Then, neglecting the second term in Eq. 
(since P2 = for |r + < a), one obtains 
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n^^^^{r,r + xz) = -9{a - b)9{z + X + w)e{w - z). (29) 

Substituting Eq. (^) into Eq. (^) we find that in the hmit A oo this equation coincides 
with our result, Eqs. (|23|) , ([26|) (assuming translational invariance for the density matrices 
Pi, P2)- We thus obtain that the GR formula reproduces the asymptotic limit, Eq. @. 
Since our result for the FSI correction factor is the exact one only for x — > oo, the same 
is true for the GR formula. Also the use of the product form for the density matrix p^, 
which takes into account only the correlation between the struck and spectator particles 
and neglects others [^, would reproduce Eq. (|2^) and therefore the GR formula (pT]). In 
this connection it is interesting to note that the same product form [^] for pA, is essential in 
an alternative derivation of the GR result for regular interactions using a modified eikonal 
formula 0, [ pJ]] . 

Another treatment of FSI in the response function, hard-core perturbation theory 
(HCPT) has been proposed by Silver [jl2|. One readily shows that in the asymptotic limit 



HCPT and GR lead to different results |]T2l. Indeed, HPCT response can be written in the 
form of Eq. (|27|) , where 

^{HCPT) ^ _ ^ (30) 

(which was used together with P2/ = d{r — a) in the comparison with GR), is different 
from il^'^^^ as given by Eq. (|29|) . It thus follows from the above consideration that HCPT 
does not reproduce the asymptotic limit for x ^ 00. 

Next, we discuss the application of our results to the response of Fermi systems inter- 



acting through a hard-core. Weinstein and Negele [|1^] studied numerically such systems. 
They considered the asymptotic limit for the response function and found considerable de- 
viations from the lA at |y| > kp, where kp is the Fermi momentum. We now derive from 
our approach a simple analytical expression for the asymptotic response of a Fermi gas with 
hard-core interaction, at small |y|. 

Consider the limit of Anpa^/d ^ 1 (the hard-core volume is small with respect to the 
volume per particle). In this case the Fermi gas momentum distribution, np(p) = ^9{kp — 
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p), is not essentially modified by hard-core interactions. However, the hard-core inevitably 
generates high momentum components in n{p). Therefore the application of our approach 
will be reliable in the region |y| < kp- The corresponding one-body density matrix is 

Pif{0; x) = 3p— . (31) 

We can consequently write the Fermi gas A-particle density matrix as 

A' 

PAF = -j^PiFix). (32) 
Substituting Eq. (|32D into Eq. (^) and taking the limit A —>■ oo we find 



4>f{y) = "~"R-6 / "^^^ ^"^^ exp(iyx — pax)dx 
TT JO kpx 



= As'-^ I i(y + -p)' - 4] 1- "^"'^M - (33) 

ATckj, ^' j + iap-kF) 2TTkl, ^ ' 

Notice that in the limit o = na^ —>■ 0, Eq. ( p3D gives the well known lA result for the Fermi 
gas response, 4>^f^{y) = 3{kp — j'^)6{kF — |y|)/4A;|,. 

In Fig. 2 we plot the resulting Fermi gas response, Eq. (^), for different values of 
the hard-core radius. We take nuclear matter as an example with kp = 1-33 fm~^ and 
p = 2kp/37T^. The solid line shows the lA result (a = 0), and the dot-dashed and dashed lines 
correspond to a=0.5 fm and a=0.8 fm, respectively. The ratio of the hard-core volume to 
1/p is 0.084 for a=0.5 fm, and 0.34 for a=0.8 fm. Therefore the use of the Fermi momentum 
distribution in Eq. ^is still approximately correct. We find that the response considerably 
deviates from the lA. One also notices the disappearance of the singularity in the response 
at |y| = kp due to hard-core FSI. 

It is also interesting to note the sizable effect from the hard-core FSI in the region of the 
quasi-elastic peak. For a=0.5 fm the response is reduced by about 10%. Although such a 
reduction has been the result of a pure hard-core FSI, a similar effect of FSI would be also 
expected for a strong repulsive interaction, provided the potential is much larger than the 
energy of the recoil particle for distances smaller than a. 
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In summary, we have investigated the effect of hard-core interactions on the response 
function obtained from inclusive scattering, around the quasielastic peak. Relevant models 
of the density matrix describing systems of particles interacting through a hard-core were 
constructed. By using them the response was shown to be modified and smeared as compared 
with the otherwise expected lA form in the asymptotic limit, both for bosons and fermions. 
New analytical formulae were obtained for the response of these systems reproducing quite 
generally the asymptotic behaviour of the response and that satisfy the asymptotic sum 
rules. 
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Figure Captions. 

Fig 1. Hard-core shadow configuration in space (r^) for tlie spectator particle, produced 
by tlie Green's function in Eq. (Q) (liatclied region). Square-hatclied circles show the regions 
of the hard-cores' interception. 

Fig. 2. Asymptotic response for a Fermi-gas in the impulse approximation (sohd), and 
with hard-core final state interactions corrections with core radius a=0.5 fm (dot-dashes), 
and a=0.8 fm (dashes). 
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